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We study the properties of the low-lying states at the critical point of the phase 
transition from U(3) to 0(4) symmetry in the U(4) vibron model in detail. By analyzing 

OO ! 

Q , the general characteristics and comparing the calculated results of the energy spectra 

and the El, E2 transition rates in E(3) symmetry, in potential model and the 
^ I finite boson number case in boson space, we find that the results in the potential 

demonstrates the characteristic of the classical limit at the critical point well and 
the E(3) symmetry over-predict the energy levels and under-predict the El and E2 
transition rates of the states at the critical point. However, the E(3) symmetry may 
describe part of the properties of the system with boson number around 10 to 20. 
We also confirm that the ^'^C+^'^C system is an empirical evidence of the state at the 
critical point of the phase transition in the U(4) model when concerning the energies 
of the low- lying resonant states. 
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Quantum phase transitions in mesoscopic system (system with a finite number of par- 



ticles N), such as atomic nuclei [iiisiEiiaiaEiiiiiHiEiiininiiniiisiEiiiaiisiEi 

[la [13 EQl EH |22l |23], molecules [Ml ESI EE] , atomic clusters [Z^ and finite polymers have 



recently been attracting a lot of interests. The transition in these systems are among dif- 
ferent shapes, geometric configurations, and modes of collective motions. For nuclei, it has 
been well known that there exists vibrational, 7-soft rotational, axially rotational, and other 
collective modes. The interacting bosons model (IBM, the simplest one is the U(6) model 
including s- and d-bosons) [6] has been shown to be successful in studying the properties 
of the low-lying collective states of even-even nuclei and the shape phase transitions (see 



for example, Refs. [IlElSSl[3[6lE[8l[am[l2lE3l[Tll6l[I7ll8lE0lEIlE2lE3] ). 



It has also been well known that there exhibits vibration and rotation in molecules. To 
characterize the relative motion of a dip ole- deformation in the three dimensional space and 
describe the behavior of rotational and vibrational motions of molecules, the U(4) vibron 
model has been developed [2ll EH] and apphed to two-body (or two-cluster) systems such as 
diatomic molecules [281 El], binary clusters [301 ED ESI [33|, qq mesons [Ml ES], and so on. It 
has been shown that the U(4) model involves two dynamical symmetries, namely U(3) and 
0(4), and there exists a second order phase transition from U(3) to 0(4) symmetry [361 137]. 
For the properties of the states at the critical point of the vibrational to 7-soft rotational, 
the vibrational to axially rotational phase transition of nucleus, respectively, there have had 
thorough studies in theory (see, for example, Refs. [TniEHlE9lllQlliTlll2lll3llll]) and quite 
a lot of empirical evidences have also been found (see, for instance, Refs. [ISl HE]). For the 
properties of the states around the critical point of the U(3) and 0(4) phase transition in 
the U(4) model, even though they can be described approximately with the E(3) symmetry, 
which has been discussed in some sense [37], they have not yet been investigated in detail. 
We will then study the low- lying energy spectrum and the El, E2 transition rates of the 
low-lying states at the critical point of the phase transition in the U(4) vibration model in 
this paper. 

In the U(4) model, elementary excitations are dipole p-bosons with spin and parity 
J'^ = 1^ and scalar s-bosons with J'^ = 0"*". With assumptions that the total number of 
bosons and the angular momentum of the system are conserved, there are only two dynamical 
symmetry limits, U(3) and 0(4). Accordingly, there exist two dynamical symmetry chains: 



It has been shown that the U(3) symmetry corresponds to nonrigid ro-vibrations, while the 
0(4) symmetry represents rigid ro-vibrations [28]. A general Hamiltonian of the U(4) vibron 
model with only one- and two-body interactions being taken into account can be expressed in 
terms of the linear and quadratic invariant operators (Casimir operators) of all the subgroups 



f/(4) D U{3) D 0(3), 
f/(4) D 0(4) D 0(3). 



(I) 
(11) 



(1) 
(2) 
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contained in the dynamical group chains. 

To study the property of the phase transition between U(3) symmetry and 0(4) sym- 
metry, one starts usually from the simple Hamiltonian 

H = e[{l-v)n-j^^D.D], (3) 

where e is a scale parameter and it can be taken as one for convenience without any loss of 
generality, n = J2mPmPm is the number operator of p-bosons, = (stp + pt§)(i) is the 
electric dipole operator, where s = s and pm = {—lY~'^p-m- f{N) is a linear function of 
total boson number N. t] is the control parameter. It is easy to show that such a Hamiltonian 
can be alternatively written as 

H = eil- v)C,ui3) - ^y^C'20(4) + ^y^^20(3) • (4) 

It is obvious that the system is in U(3) symmetry when rj = 0, and in 0(4) symmetry if 
1] = 1. By varying t] G [0, 1], we can realize the U(3)-0(4) phase transition. 

The classical limit corresponding to the Hamiltonian in Eq. ([3]) can be obtained by 
considering its expectation value of the coherent state 



\N- 1) = (iV!)-i/2[(i _ t* . ty/^s^ + t . pt]^|0) , (5) 

where t is a complex three-dimensional vector, and its complex conjugate is denoted by 
t*. Then the classical Hamiltonian can be given as Hd = {N;ti \ H \ N;t). One can 
introduce canonical position and momentum variables r and q (they can be the ones in 
three dimension) by the transformation [24J 



t = (r + zq)/V2, t* = (r-iq)/V2. (6) 

So the classical potential is just the value of Hci{q,r) with q = 0, where |r| = r and |q| = q 
i.e., 

V{r) = H4q = 0,r). (7) 

In the case of taking f{N) in Eq. ([3]) to be 3A^, the energy potential corresponding to the 
Hamiltonian in Eq. ([3]) can be concretely expressed as 

V(r) = iV[(l-r/)^- V(2-r2)]. (8) 

From Eq. (jS]), we could easily recognize that the classical potential for U(3) symmetry 
(with 1] = 0) is V[/(3)(r) = ^r^, while the classical potential for 0(4) symmetry (with i] = 1) 
is Vo(4)(r) = — yr^(2 — r^). By analyzing the stability of the system with potential in Eq. ([8]) 
we could find that the critical point of the quantum transition corresponds to the control 
parameter r]c = f . It is obvious that, as r/ < |, Vmin{r) = 0; if r/ > |, Vmin{r) = 
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Furthermore, at the critical point with control parameter r] = rjc = the energy potential 
can be explicitly written as 

VcAr) = , (9) 
and the Hamiltonian corresponding to Eq. ^ could be simplified to a Schrodinger equation 



"1 A 



$ = E$. (10) 



Looking through the characteristic of the potential at the critical point in Eq. ([9]), one can 
learn that it is quite fiat around its bottom. To study the property of the states at the critical 
point of the phase transition in the U(4) model, at first one can go along the way taken in 
Ref. [38] and approximate the potential around the critical point to be a three-dimensional 
infinite square well 

nr) = l°' (11) 
I oo , r > rw ■ 

It is apparent that such a potential deviates from that of the E(5) symmetry of the transition 
from U(5) to 0(6) symmetry in the IBM [38] only in dimension. Then the states generated 
from this potential can de denoted as the ones with the E(3) symmetry. 

In such a situation, the Schrodinger equation is exactly solvable and the solution could 
be expressed as some Bessel function. The excitation energy can be given as 

K,L = ^<L, (12) 

with k"^ ]^ = where yn,L is the nth zero point of the Bessel function Jl+i/2{z), B is a 
constant. We give in Fig. [1] the obtained energy spectrum of some of the low- lying states 
and in Table [1] the values of the excitation energies of the states, where the energy of the 
ground state is set to zero and all energies are normalized to the energy of the first excited 
state (with = li). 

Then we solve the Eq. (fTOl) numerically to discuss more practically the properties of the 
low-lying states at the critical point of the phase transition in the U(4) model. The obtained 
energy spectrum of the low-lying states is illustrated in Fig. [2] and the concrete values of 
some states' energies are listed in Table [H To show the variation feature of the energy of 
low-lying states explicitly and make it easy to compare with experimental data, we display 
some of the energy ratios between some states in Fig. [3l Comparing the result in the E(3) 
symmetry with that obtained by solving the Schrodinger equation in potential of form, 
one can notice that, the energy of a state in the E(3) symmetry is higher than that given by 
the potential for the state with the same quantum number. 

To investigate the properties of the states at the critical point more comprehensively, 
we also solved the engin-equation with Hamiltonian in Eq. ([3]) with t] = j directly in the 
cases of various boson numbers. The obtained results of the variation feature of the ratios 
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Figure 1: Calculated energy spectrum of the low-lying states in E(3) symmetry (with a 
potential of infinite square well) around the critical point of the phase transition U(3)-0(4) 
in the U(4) vibron model and the corresponding El (panel (a)), E2 (panel (b)) transition 
rates (the numbers close to the arrow are the value of the rate with normalization -B(E1; 1^ — > 
0^) = 100, B{E2] 2t Of) = 100, respectively). 



of some low-lying states against the boson number and the comparison with those in the 
E(3) symmetry and potential are illustrated in Fig. HI From Fig. HI one can learn that, 
for the yrast states, the results given in the E(3) symmetry is quite close to that in the case 
of fewer bosons (for instance, around ten), but that determined by the potential is the 
asymptotic limit of the one in large boson number limit. It indicates that the one with the 
potential describes the behavior at the classical limit (or infinite boson number limit) 
more appropriately, and the E(3) symmetry over-predicts the energy levels. Such a feature 
is consistent with those of the states around the critical point of the phase transition from 
U(5) to 0(6) in the U(6) model of IBM ^ Sg. 
With transition operators 

f (El) = -P-Eo = ei,e//rcos0, (13) 
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Figure 2: Calculated energy spectrum of the low-lying states at the critical point of the phase 
transition U(3)-0(4) (with a potential of r^) in the U(4) vibron model and the corresponding 
El (panel (a)), E2 (panel (b)) transition rates (the numbers marked close to the arrow are 
the value of the rate with normalization B{E1; 1^ Of) = 100, B{E2;2t Ot) = 100, 
respectively). 



1 f)J^ 

f (E2) = --T.Q4^), = e,,effr\3cos'e- 1) , (14) 

where Qij = q{3xiXj — r'^Sij), we can determine the El and E2 transition rates B{Ek, Li 
Lf) = 2^q:j-|(n/ L/||T(Ek)||nj Lj) p. The obtained results in case of the E(3) symmetry (i.e., 
with the potential of infinite square well) and those in potential of are displayed in Fig. [H 
Fig. [21 respectively. We have also calculated the El and E2 transition rates with operator 
T{El)q = ei,e//(s^p + p^sYq, T{E2)^ = e2,e//(p^p)^/^, respectively, for the systems with 
different boson numbers. The calculated results of changing behavior of the ratios of some 
transition rates with respect to the boson number and the comparison with those in E(3) 
symmetry and potential are manifested in Fig. [51 Looking through Figs. [1], [2] and [5l one 
can notice that the E(3) symmetry under-predicts the El and E2 transitions rates between 
the low- lying states at the critical point with the potential in the U(4) model. Moreover, 
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Table 1: Excitation energies of the low- lying states at the critical point of the phase transition 
from U(3) to 0(4) in the U(4) vibron model. 
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the E(3) symmetry can only describe the ratios of some transition rates of the system with 
boson number N around 10 to 20 with an exception of f(f 2-21^00 ' "^^ich deviates from both 
the result in E(3) symmetry and that in the potential in small N case but in large-A^ 
limit obviously approaches to the one in the potential. Therefore, the potential depicts 
excellently all the transition rate ratios as those energy ratios of the system in classical limit. 
The analyses of energy levels together and transition rates indicate that the results in the 
E(3) symmetry shows a nature finite N correction on those in the potential for the low- 
lying spectrum since the bottom of the potential in E(3) symmetry is flatter than that in the 
potential. As for the high-lying spectrum, this conclusion may be changed, which needs 
to be further tested. 

It has been shown that the U(4) model can successfully describe nuclear molecules [171 
HE]. A nuclear vibron model for nuclear molecules consisting of two clusters holds generally 
a dynamical symmetry Gci ® Gc2 ® Ufi{4), where the internal structure of the ith cluster 
is described by Gci which may be, for example, the U(6) IBM or SU(3) shell model, and 
the relative motion between the clusters is described by the Uij(4) vibron model. In some 
cases, only the U(4) vibron model itself is sufficient to describe the rotational and vibrational 
excitations in nuclear molecules, where the internal structure of each cluster does not play 
an essential role in the low-lying levels, such as, the narrow resonances in the ^^C+^^C 
system [HH [32]. In the early time, the 0(4) limit of U(4) vibron model was proposed to 
describe the resonant energy of ^^C+^^C system |31j, while the analysis in Ref. [32] indicates 
that the U(3) limit may be more preferred when fitting the energies of the low-lying resonant 
states, and Ref. [H] shows that the experimental data are in fact between those with U(3) 
and 0(4) symmetry, respectively. Recently some of us proposed that the ^^C+^^C system 
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Figure 3: Calculated ratios of the energies of some low-lying states '^(2")^^(o") -E(Oi) = 

and n = 1 panel (a), = 2 panel (b)) and some of the experimental data of ^^C+^^C system. 



may be described by the critical symmetry, E(3), in the U(3)-0(4) phase transition P?] . 
Since the potential at the critical point of the U(3)-0(4) phase transition is in fact in the 
form of but not that in the E(3) symmetry, to show the practical possibility of such a 
system being that at the critical point of the phase transition in the U(4) model, we should 
re-analyze that more cautiously. Then we display the experimental data of some ratios of 
the low-lying resonant energies of the ^^C+^^C system as well as the corresponding results 
with U(3), 0(4), E(3) symmetry and r^-potential in Table [2] and Fig. [31 One can find 
from Table [2] and Fig. [3] that the experimental data agree globally better with the results 
in r^-potential than with those in the U(3), 0(4) and E(3) symmetries. It indicates that, 
concerning energies of the low-lying resonant states, the ^^C+^^C system is an empirical 
evidence of the states at the critical point of the U(3)-0(4) phase transition in the U(4) 
model, and the E(3) symmetry can describe that approximately. 

In summary, we have calculated the energy spectra and the transition rates of not only 
El but also E2 transitions of the low-lying states at the critical point of the phase transition 
from U(3) to 0(4) in the U(4) vibron model in coordinate space with both the E(3) symmetry 
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Figure 4: Calculated variation behavior of the energy ratios of some low-lying states with 
respect to the boson number and those in the E(3) symmetry and in the potential. 
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Figure 5: Calculated variation behavior of the ratios of some El, E2 transition rates between 
low-lying states with respect to the boson number and those in the E(3) symmetry and in 
the potential. 



(or with potential in the form of infinite well) and the potential as well as in boson space 
with finite A^. Our calculation shows that the large- limit of the U(4) vibron model at the 
critical point of phase transition can be represented by a r"^ potential model excellently, but 
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Table 2: Experimental data of the energy ratios of some low-lying resonant states of ^^C+^^C 
and the corresponding values in the r^-potential, E(3) symmetry, U(3) symmetry and 0(4) 
symmetry (the experiment data are taken from Ref. [IE])- 





Expt. 


E(3) 


^4 


U(3) 


0(4) 




2.22 


2.45 


2.20 


2 


10 
3 




5.72 


4.31 


3.55 


3 


7 


-^02 /-^2i 


0.31 


1.27 


1.05 


1 


2N 
3 


-^22 7-^21 


2.50 


3.12 


2.39 


2 


3 ^ 




4.48 


5.45 


3.79 


3 


2N 1 10 
3 3 




7.55 


8.22 


5.31 


4 


^ + 7 
3 ^ ' 




1.86 


3.38 


2.47 


2 


4iV-4 
3 


E^^/Eq^ 


7.16 


1.93 


1.97 


2 


5 

N 


Eq^/Eq^ 


18.45 


3.39 


3.19 


3 


21 
2N 


Ea^/Eq^ 


1.19 


0.72 


0.89 


1 


5 

2iV-2 


E&^/Eq^ 


3.08 


1.27 


1.43 


1.5 


21 

4Ar-4 



not close to those of an infinite well in the i^^(3) symmetry. Generally, The E(3) symmetry 
over-predicts the energy levels and under-predicts the El and E2 transition rates of the 
states at the critical point of the phase transition in the U(4) model in large-A^ limit, but 
predicts a critical-point spectrum that is qualitatively similar to the U(4) vibron model for 
small values of in the most cases. The calculated results also express that the ratios 
of quantities at the critical point quickly approach to a constant with the increasing of 
boson number N when > 30, which further confirms the N scaling behavior of quantities 
at the critical point shown in Ref. [37], where the results indicate that the ratios of both 
energies and El transition at the critical point are approximately invariant as N varies. 
Comparing the theoretical results with the experimental data of some systems, we find that 
the ^^C+^^C system is an empirical evidence of the state at the critical point of the U(3)- 
0(4) phase transition in the U(4) model and the E(3) symmetry can describe that roughly 
when concerning the energies of the low-lying resonant states. To conform it much more 
solidly, one needs the data of electromagnetic transition rates. 
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